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OitMarazfN.: 


Abstr  art 


We  develop  new  confidence  interval  estimators  for  the  underlying  mean 
of  a  stationary  simulation  process.  These  estimators  can  be  viewed  as 
generalizations  of  Schruben's  so-called  standardized  time  series  area 
confidence  interval  estimators.  Various  properties  of  the  new 
estimators  are  given. 


In  thi*.  note,  we  present  new  confidence  interval  estimators  for 


the  underlying  mean  u  of  a  stationary  simulation  process.  These 


estimators  can  be  viewed  as  generalizations  of  the  so-called  area 


estimators  given  in  Schruben  (1983)  and  Goldsman  (1984). 


Consider  the  stochastic  process  X,,...,X  .  Define 

I  m 


VJ~j  Ek=l  Xk  ’ 


=  lim  mVar(X  )  ,  and 

m-*®  m 


T  (t)  = 
m 


|mt | (X  -X.  .  .  ) 
L  J  m  |_mt  J 


,  °< t< 1 , 


where  I.  ]  is  the  greatest  integer  function.  tT  (t)  ,  0<t<l>  is 

“  *  m  —  — 


called  the  standardized  time  series.  Schruben  (1983)  proves  that  if 


X^,...,X  is  a  stationary,  e-mixing,  finite  variance  sequence  of 


random  variables  (satisfying  one  other  technical  condition),  then  as 


m-«®,  T  (t)  converges  in  distribution  to  a  standard  Bronnian  bridge 


process,  fB^,  0<t<_l>.  Also,  the  standardized  time  series  is 


asymptotically  independent  of  mX 


Remark;  It  is  well  known  that  B.  ~  Nor (O, t ( 1-t) )  and  Cov(B.  ,B  )  = 

t  tj  t2 


mi  n ( t  ^ , t^ )  ( 1 —max (t^,t2)). 


Define  A  s  ^  ’  the  ck  s  be*n9  pre-specif i ed. 


For  large  m,  A  :=•  (is  approximately  distributed  as)  <rNor(0,V),  where 


V  =  1=  J"  Cov(c  .B.  .  ,c.B.  ) 

2  =  l  wk=l  j  j/m’  k  k/m 


^W5«<S<lwSWS' 


mm 


It  will  be  computationally  convenient  to  approximate  V  by  letting 


j 


I 

I 


* 

I 


ms  and  k  -  mt  (so  that  djdk  ®  m*"dsdt).  Thus, 

V  •=  cjck  m*n  < »k/n*>  1 1-max  ( j/m,k/m)  3 


c  c  .  min (s,t> C 1— max (s, t ) 3 
ms  mt  ’ 


.  c  c  .  s(l-t)  dsdt 
0  ms  mt 


dsdt 


(1) 


2  2  2 

So  A  / V  ■*  a  %  (1),  and  this  is  called  the  weighted  area 

2 

estimator  -for  the  variance  <r  . 

Suppose  now  that  we  work  with  the  process  X_,...,X  ,  where  n=bm, 

l  n 

and  that  this  series  satisfies  Schruben's  conditions.  Divide  the 

process  into  b  contiguous  batches ,  each  of  size  m;  i.e.,  X..  ,  , 

X(i-nm+2’  •••»  Xim  comprise  batch  i,  i=l,...,b.  Each  individual 

batch  can  be  standardized:  For  isl,...,b  and  j:=l,...,m,  let 

y  _  cj  „  (average  of  the  first  j  X's 

i»j  -  j  £k=l  (i— l)m+k  from  the  i-th  batch). 


(grand  mean) , 


T. 

l 


(t) 


,m 


I  mt  I  (X .  -X.  ,  ..) 

L  J  i  ,m  l , [mt J 


eVm 


0<t£l ,  and 


A. 

l 


c  r-m 
m  wk=l 


c. T.  (k/m)  . 
k  l  ,m 


For  large  enough  m,  each  of  the  standardized  time  series  Cthe 
T  (t) *s3  is  approximately  distributed  as  a  Brownian  bridge;  so  A. 

l  ,m  i 

<rNor(0,V),  i=l,...,b.  Further,  for  large  m,  we  can  treat  the  batches 


as  if  they  were  (approximately)  independent.  This  yields: 


An  immediate  consequence  of  Theorem  21.1  of  Billingsley  (19&U)  is  that 


x-u 

Z  =  — - -  :=•  Nor  (0,1)  . 

n  a /Vn 


Then  the  asymptotic  independence  result  above  gives: 


n 


[L!tr2 

L  <r2bV  • 


Nor (Qt 1 ) 

[^]‘/2 


-  t  (b) 


We  -finally  obtain  confidence  interval  estimators  for  u: 


r  .  1/2 

Pr  {  w  £  Xn  !  tb,l-a/2^  nbV  £i  =  l  Ai  ^  } 


1-a  , 


(2) 


where  t^  ^  is  the  upper — v  quantile  of  the  t(b)  distribution. 

We  consider  various  choices  for  the  weights  (the  c^'s).  These 
choices  and  their  resulting  V's  Cfrom  the  integral  approximation  (1)1 
are  summarized  in  Table  1. 


Remarks: 

(i)  After  choosing  the  weighting  sequence  <c^>,  the  associated  V  from 
Table  1  is  used  in  (2)  to  form  a  confidence  interval  estimator. 

(ii)  The  calculation  of  V  from  (1)  is  straightforward  (but  sometimes 
tedious) . 

Example:  For  choice  3  (from  Table  1), 


-A- 


v  !  r  Ji  Jo  6<I-t>  d‘  dt 


^  fl  ft  1  ^  „  fl  1  r»<l-t>  ^ 

-  Jo  Jo  TPTTt  d5  dt  "2  Jo  t -  dt 


(lii)  The  variance  estimator  which  results  -from  choice  1  <the  equal 
weighting  case)  is  asymptotically  the  same  (as  m-*»)  as  the  so-called 
area  variance  estimator  from  Schruben  (1983). 

(iv)  For  each  standardized  time  series,  CT.  (m)>,  choice  2  grants 

i ,  m 

greater  weight  for  'small '  values  of  t.  Choices  3  through  6  give 
comparatively  little  weight  to  the  middle  (t  :=•  1/2)  of  each 
standardized  time  series. 


(v)  Table  2  summarizes  an  empirical  study  involving  the  order  1 
exponential  autoregressi ve  model  Ccf.  Lewis,  (1980)3.  The  weighted 
area  estimators  are  seen  to  perform  well  for  'large'  batch  size. 


(vi )  Denote  the  random  variable  corresponding  to  the  half-length  of 
the  weighted  area  estimator  by  H.  Following  Schmeiser  (1982)  and 
Goldsman  and  Schruben  (1984)  (G-S) ,  it  is  easy  to  derive  the 

f ol lowing: 


ECH3 


<r 

Vn 


l-a/2 


(2/b) 


1/2  r((b+l)/2) 
r(b/2> 


* 


Var (H) 


-ib 
n  b , 


"  1  -  2 

■  r((b+l)/2)  i2  1 

a/2 

L  1  & 

L  r(b/2)  J  J 

and 


The  coverage  probability. 


Pr<‘W  <  "  F(tb,i-a/2>  "  F<‘tb,l-a/2) 


where 


r<.)  is  the  gamma  function  and  F(.)  is  the  c.d.f.  of  the  noncentral 


-6- 


of  correlation  amongst  batches  which  are  encountered  when  using  the 
area  estimator. 
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Table  2 


(•'erlormance  of  weighted  area  confidence  interval  estimators  for  the 
mean  of  an  CAR  ( 1 )  process  with  coefficient  />  =  0.2  and  exponential 
<mean=l)  noise  based  on  100  independent  runs  of  2560  observations 
each.  [Choices  of  weights  are  summarized  in  Table  13. 


weighting  | 

choice  |1  2  3  4<e=0)  5<e=0) 

I 
I 

Confidence  interval  achieved  coverage  <907.  desired) 


b  m 


1 

2560 

1  .88 

.92 

.90 

.89 

.90 

2 

1280 

1  .93 

.96 

.94 

.93 

.95 

5 

512 

1  .95 

.94 

.96 

.95 

.96 

10 

256 

1  .93 

.93 

.94 

.95 

.96 

20 

128 

I  .93 

.94 

.91 

.91 

.91 

Average 

confidence 

interval 

half 

-lenath  (  x 

10000) 

1 

2560 

1 

1 

1  1193 

1266 

1184 

1199 

1202 

2 

1280 

1  572 

612 

598 

613 

617 

5 

512 

1  480 

481 

486 

492 

498 

10 

256 

I  440 

441 

436 

444 

446 

20 

128 

I  426 

423 

415 

427 

426 

U'-' 


Sample  standard  deviation  of  half-lengths  (  x  10000) 


1 

2560 

I  900 

923 

876 

932 

931 

2 

1280 

I  297 

300 

298 

294 

298 

5 

512 

I  153 

149 

150 

150 

149 

io 

256 

l  103 

104 

97.5 

95.6 

93 

20 

128 

|  52 . 6 

50.7 

57.9 

59.5 
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